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The Fundamental Plane of Spiral Galaxies: 
Search from Observational Data 
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Abstract 

The fundamental plane of spiral galaxies was searched from observational data, which can be approxi- 
mately represented in observational parameters by i oc V^R, where L , R, and V are the luminosity, the 
linear size of galactic disk, and the rotation velocity. This plane exists at all optical bands in our samples of 
more than 500 spiral galaxies in total. It is more fundamental than the relationship between any two of the 
parameters, and can reduce the residual of the TuUy-Fisher L-V relations by about 50%. Noticed that the 
power index of V is doubled from that of R, which implies that the total mass and the mass distribution 
of galaxies plays an important role in forming the fundamental plane. Involving the third parameter of 
galactic size has a strong physical implication on galaxy formation and the dark-mass distribution. 

Key words: galaxies: fundamental parameters — galaxies: kinematics and dynamics — galaxies: 
spiral — galaxies: statistics 



1. Introduction 

The TuUy-Fisher relation is a tight correlation between 
the internal motion and the luminosity of spiral galaxies 
(TuUy, Fisher 1977). It can be expressed as 



M^alogV + 'j, 



(1) 



where M is the absolute magnitude, V the rotation veloc- 
ity, a the slope of the relation, and 7 the zero-point. This 
empirical relation has been used for estimating the dis- 
tances of galaxies, and hence for determining the Hubble 
constant (e.g. Sakai et al. 2000; TuUy, Pierce 2000). It 
also provides a critical constraint on galaxy formation (e.g. 
Dalcanton et al. 1997; Mo et al. 1998; van den Bosch 
2000; Avila-Reese, Firmani 2000; Firmani, Avila-Reese 
2000; Koda et al. 2000b). However, the detailed phys- 
ical processes concerning the origin of the TuUy-Fisher 
relation have not yet been fully understood. This rela- 
tion could be a result of the self-regulated star formation 
in the disks of different mass (see Silk 1997; Eisenstein, 
Loeb 1996), or a direct consequence of the cosmological 
evolution between the mass and the circular velocity (see 
e.g. Stcinmetz, Navarro 1999; Mo et al. 1998; Mo, Mao 
2000). 

It was believed that the residual value deviating from 
the TuUy-Fisher relation can pose strong constraints on 
scenarios of galaxy formation and evolution (see e.g. 



Eisenstein, Loeb 1996) which was systematically stud- 
ied by, for example, Willick ct al. (1997). The residuals, 
written as 

S = M-{a\ogV + -f), (2) 

are not accountable by measurement errors. With good 
data at the R and / bands, the scatter of the lumi- 
nosity may be partially attributed to the measurement 
uncertainty of 0.08 mag (Courteau 1996) or 0.04 mag 
(Giovanelli et al. 1997), which is further convolved with 
worse-determined distances. Because of the residuals, the 
estimates of the galaxy distance based on the TuUy-Fisher 
relation have a typical uncertainty of 20%. 

Some efforts have been made previously to search for 
tighter correlations among the luminosity, rotation veloc- 
ity, and disk radius for spiral galaxies. Kodaira (1989) 
found a much tighter correlation among these three pa- 
rameters, L oc VR^. He pointed out that the above rela- 
tion is valid for ellipticals as well. Very recently, Koda et 
al. (2000a) found that the /-band luminosity L, /-band 
radius /?, and the rotation velocity V of spiral galaxies 
are distributed on the so-called scaling plane in three- 
dimensional parameter's space of these quantities, which 
can be expressed as L oc {VRY'^. Because the galactic 
radius was considered (Koda et al. 2000a), the residual 
of the TuUy-Fisher relation was significantly reduced. In 
fact, Courteau and Rix (1999) also realized the impor- 
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tance of involving the disk size to the TuUy-Fisher re- 
lation, but they did not come up with such a relation. 
Willick (1999) explored the role of the disk scale length, 
and found the surface brightness dependence of the TuUy- 
Fisher relation. 

It is intriguing to know whether and how the third pa- 
rameter, the galactic size, really plays some role in galaxy 
formation, and how the mass and luminosity are physi- 
cally related to this parameter. Besides strong evidence 
found by Koda et al. (2000a), the linear (optical) galactic 
diameter has a very tight correlation to the H i mass in 
later galaxies (Broeils, Rhee 1997). Koda et al. (2000b) 
have already done a set of simulations. They found that 
the galactic size was involved as a key parameter in such 
a process that the galactic mass and angular momentum 
were controlled during galaxy formation. Very recently, 
Shen et al. (2001) considered the variations of the model 
parameters presented in the theoretical work of Mo et al. 

(1998). They also found a theoretical fundamental plane, 
Locy2.6^o.5^ 

It is therefore very necessary to search the fundamental 
plane from observational data with the galactic radius, 
R, as the third parameter, and see how it matches the 
results from simulations of galaxy formation. The primary 
procedure we take here is to see how the galactic size can 
help to reduce the residual of the TuUy-Fisher relation. 
The best-fitting plane, expressed as 



M ^ alogV + l3logR + J, 



(3) 



was searched from observational data with three variables, 
a, (3, and 7. Koda et al. (2000a) obtained a = (3 = —3.25, 
while Kodaira (1989) gave a == -2.5 and /3 = -5.0. We 
found that the best-fitting plane with fi/a around 0.5 for 
different data-sets, namely, the /-band data of Han (1992) 
and Palunas and Williams (2000), the i?-band data of 
Courteau (1996, 1997), and the BVRIH-h&nd data of 
calibration galaxies in Sakai et al. (2000) and Macri et al. 
(2000). The plane exists in different wavebands, and can 
reduce the residual of the TuUy-Fisher relation by about 
50%. This implies that about 50% of the scattering in 
the observed TuUy-Fisher relations is dominated by the 
intrinsic properties of galaxies, and the other 50% by mea- 
surement uncertainties. We believe that the plane is fun- 
damental for spiral galaxies, similar to that for elliptical 
galaxies (Djorgovski, Davis 1987; Dressier et al. 1987). 
It should provide useful constraints on models of galaxy 
formation (see Shen et al. 2001; Koda et al. 2000b). 

2. The Fundamental Plane versus the ScaUng 
Plane 

Han (1992) presented carefully-calibrated /-band data 
for galaxies in clusters, including the total magnitude, 
M)°* and M]^-^ (in units of mag), and the face-on /-band 
isophotal radius, A^^ 5 (in units of arcsec). He also gave 
the width of the 21 cm Hi line, W20 (in units of km s~^) 
corrected for the redshift and inclination. The scaling 
plane, L oc {VRf-^ (see figure 1 of Koda et al. 2000a), 
was found among these three observational parameters. 



Therefore, the TuUy-Fisher relation between A// and W20 
is an oblique projection of the plane. This also explains 
the Freeman relation, L oc R^ (Freeman 1970), and the 
radius-velocity relation (TuUy, Fisher 1977). 

We used the same /-band data-set to fit equation (3). 
In practice, we define k = f3/a and fit equation (3) by 
varying a and k. Obviously, equation (3) becomes the 
original TuUy-Fisher relation when k = 0.0, or it goes to 
the scaling plane of Koda et al. (2000a) when k = 1.0. 
After excluding those galaxies (1) whose recession veloc- 
ities deviate by more than 1000 km s~^ from the mean 
velocity of a cluster, (2) without Hi measurements and 
(3) which may not be a cluster member, we finally ob- 
tained 160 galaxies in the sample. Throughout this pa- 
per Hq = 75 km s^^Mpc^^ is assumed. We first used 
the /-band total absolute magnitude data, M = M)°* , to- 
gether with V = W20 and R = R123.5 (in kpc), and ob- 
tained the minimum x^= 252.2 at (a,/3,7) = (-4.05±0.21, 
-3.00 ±0.14, -8.01 ±0.43), see figure 1. This minimum 
X^ value is much smaller than the minimum x^ — 697.6 
from the TuUy-Fisher relation (fc = 0). It can also be seen 
in figure 1 that the data scatter is much smaller when 
the galactic radius is involved. The best-fitting value of 
fc = /3/Q;is0.74±0.06, i.e. it is at neither A: = 0.0 nor fc = 1.0 
(figure 1). All fitting results, including those found below, 
are listed in table 1 . The last two columns are the x^ val- 
ues and the scattering of luminosity S, both for the original 
TuUy-Fisher relation and for the best fitting plane. 

To account for light within the linear radius, the ab- 
solute magnitude, M = Mf^-^, should probably be used. 
We then obtained the best-fitting parameters: (a,/3,7) = 
(-4.75±0.21, -2.85±0.14, -6.21±0.43) at fc = 0.60±0.04 
and x^ = 359.4. 

3. The Fundamental Plane in Other Bands 

As a matter of fact, not many publications about the 
surface photometry of spiral galaxies contain information 
about galactic sizes. The database of TuUy-Fisher cali- 
bration galaxies (Macri et al. 2000; Sakai et al. 2000), 
the database of Courteau (1996, 1997), and the very re- 
cently published /-band data by Palunas and Williams 
(2000) were used to search for the fundamental plane. 

3.1. The Hubble Calibration Galaxies 

Sakai et al. (2000) pubhshed the Cepheid distances 
and carefully-corrected absolute magnitudes at BVRIH- 
bands of 21 calibration galaxies for the Hubble Space 
Telescope (HST) key project for the distance scale. They 
also presented the well measured and corrected Hi 20% 
line widths. Macri et al. (2000) further published the 
brightness profiles at the BVRI-haiiAs and isophotal radii 
Rb25 and /?723.5- We measured the the isophotal radii 
{Rr23, Rr25, Rv2z, Rv25, Rb23, and Ri2b) from the 
brightness profiles. Finally, there are 17 galaxies with all 
parameters which we needed. 

We fitted equation (3) to these observational parame- 
ters in each band. As expected, involving the radius can 
greatly reduce the scattering around the TuUy-Fisher re- 
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Fig. 1. Comparison of the TuUy— Fisher relation (the left side) and the best fitting plane (the right side). Data of 160 galaxies were 
taken from Han (1992). The typical rms errors of the absolute magnitudes, together with the displacement from the uncertainties 
of W20 and -R23.5, were taken as 0.2 mag. Galaxies in different clusters were plotted by different symbols. The data distribution in 
3-D space and the variation of fitting x^ with k were also plotted. 



lation (figure 2). The k values are around 0.5 (see lines 
with HST- in table 1). Notice that in these HST data- 
sets, those with R123.5 and Rb25 produced the smallest 
residuals. 

3.2. I -Band Data of Palunas and Williams (2000) 

Very recently, Palunas and Williams (2000) conducted 
/-band photometry and Ha measurements of 74 galaxies. 
They presented the corrected magnitudes, M)°* and radii, 
R123.5, at the 23.5 mag arcsec"^ isophote. They showed 
that for Freeman type-I galaxies, i.e., galaxies with an 
exponential disk, the disk scale length, r^, is very well- 
correlated with the disk size of R123.5, while this is not the 
case for Freeman type-II galaxies with a large fraction of 
disk which does not follow the exponential disk distribu- 
tion (see their figure 2). For the TuUy-Fisher relation, 
they measured the rotation speed from the weighted aver- 
age of the rotation curve points, and defined the velocity 
width as twice the speed. We corrected these widths for 



the inclination angle and redshift effect. 

As can be seen in figure 3, we found that the fitting 
with R123.5 can largely diminish the scattering around the 
TuUy-Fisher relation. See table 1 (see lines with PWOO) 
for the fitting results. When the disk-scale length, r^, 
was used, the residual could be maximally reduced when 
fc = 0.35. 

3.3. R-Band Data of Courteau (1996, 1997) 

Courteau (1996, 1997) published the kinematic (from 
optical spectroscopic observations) and photomatic (at R- 
band) data of 304 late- type spiral galaxies. Besides the 
isophotal radii, Rr23 and Rr25, he also obtained a scale 
length for exponential disks, rj, to express the disk size. 
Courteau (1997) found that the peak rotational velocity, 
V2.2, measured at i? = 2.15rd matches the 21 cm Hi line- 
widths best, and yields the smallest TuUy-Fisher residu- 
als. Therefore, we use the velocity, V2.2, for our discussions 
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Table 1. Search for the fundamental plane from observational data. 



Sources 


Input data 
L V R 


Number 
of objects 




Fitting results 








(5tf/5 




a 


P 


7 


k 


X^f/x^ 


(mag) 


Han92-I-1 


Mr 


W20 


R123.5 


160 


-4.05 ±0.21 


-3.00 ±0.14 


-8.01 ±0.43 


0.74 ±0.06 


697.6/252.2 


0.42/0.25 


Han92-I-2 


Mf-5 


W20 


R123.5 


160 


-4.75 ±0.21 


-2.85 ±0.14 


-6.21 ±0.43 


0.60 ±0.04 


760.0/359.4 


0.44/0.30 


HST-I-1 


j^tot 


VK20 


R123.5 


17 


-5.67 ±0.38 


-2.88 ±0.29 


-3.92 ±0.80 


0.51 ±0.06 


124.8/ 


17.8 


0.32/0.15 


HST-I-2 


Mf* 


W2Q 


R125 


17 


-6.84 ±0.34 


-2.04 ±0.25 


-1.63±0.74 


0.30 ±0.04 


124.8/ 


56.8 


0.32/0.24 


HST-R-1 


MT 


W20 


Rr23 


17 


-5.30 ±0.38 


-2.26 ±0.25 


-5.37±0.81 


0.43 ±0.06 


113.7/ 


29.2 


0.32/0.17 


HST-R-2 


MT 


W20 


Rr25 


17 


-6.19±0.32 


-2.10 ±0.25 


-2.75 ±0.68 


0.34 ±0.04 


113.7/ 


42.2 


0.32/0.21 


HST-V-1 


MJ,°* 


W20 


Rv23 


17 


-4.16±0.41 


-2.54 ±0.23 


-7.86±0.91 


0.61 ±0.08 


144.7/ 


23.0 


0.35/0.16 


HST-V-2 


M'y°' 


W20 


Rv25 


17 


-5.36 ±0.34 


-2.59 ±0.25 


-4.04 ±0.74 


0.48 ±0.06 


144.7/ 


32.9 


0.48/0.19 


HST-B-1 


MT 


W20 


Rb23 


17 


-3.85 ±0.44 


-2.30 ±0.26 


-8.62 ±0.99 


0.60 ±0.10 


124.5/ 


47.0 


0.41/0.25 


HST-B-2 


MT 


W20 


Rb25 


17 


-4.19±0.39 


-3.41 ±0.33 


-5.69 ±0.85 


0.81 ±0.10 


124.5/ 


18.9 


0.41/0.18 


PWOO-I-1 


M)°* 


2{V) 


R123.5 


74 


-2.56 ±0.26 


-3.96 ±0.24 


-11.1 ±0.51 


1.55 ±0.18 


371.8/ 


88.5 


0.51/0.22 


PWOO-I-2 


Mf* 


2{V) 


Td 


74 


-5.16±0.19 


-1.79±0.16 


-7.88 ±0.45 


0.35 ±0.03 


371.8/252.9 


0.51/0.38 


C96-R-1 


Mr 


V2.2 


2.15rd 


285 


-4.11 ±0.19 


-2.23 ±0.12 


-8.45 ±0.41 


0.54 ±0,04 


703.5/359.4 


0.47/0.34 


C96-R-2 


MT 


V2.2 


Rr23 


285 


-2.25 ±0.23 


-3.64 ±0.16 


-11.4±0.46 


1.62 ±0.18 


703.5/] 


183.3 


0.47/0.24 




2.4 



2.6 2.8 2.8 3 3.2 3.4 



;(■== 144.7 




17 Galaxies 



_J 1 I 1 I 1 L_ 



2.4 2.6 2.8 2.8 3 3.2 3.4 

log 1^20 logr2„+0.61*logi?^3 



C\2 - 

I 



!^ C\3 



= 113.7 




17 Galaxies 



2.4 2.6 2.8 2.8 3 3.2 

logW^ logB^20+0-43*logi?B23 




17 Galaxies 



2.6 2.8 3.2 3.4 3.6 3.8 4 

log 1^30 logW^^+OmHogRg^^ 



Fig. 2. Scattering from the TuUy-Fisher relation (left) and the best-fitting plane (right) for the Hubble calibration galaxies at the 
IRVB bands. Data were taken from Sakai ct al. (2000) and Maeri et al. (2000). 
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Fig. 3. Comparison of the TuUy— Fisher relation (the left side) and the best-fitting plane (the right side). Data of 74 galaxies were 
taken from Palunas and Williams (2000). The typical rms errors of the absolute magnitudes were taken as 0.2 mag. The data 
distribution in 3-D space was also plotted. 
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Fig. 4. Comparison of the TuUy-Fishor relation (the left side) and the best plane (the right side). Data of 285 galaxies were taken 
from Courteau (1996,1997). The typical rms errors of the absolute magnitudes were taken as 0.3 mag. The variation in the fitting 
residual x^ with the ratio k = /3/a is plotted on the right side. 
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below. We took these parameters^, together with the cor- 
rected total absolute magnitude, M^^, from his combined 
data set (Courteau 1999, astro-ph/9903297). 

As listed in table 1 (lines with C96-), using the scale- 
length to express the disk size gives the minimum at 
k ~ 0.5, but the residual is reduced only by a factor of 
2 compared with the TuUy-Fisher relation. See figure 4 
for illustrations. In contrast, using the isophotal sizes 
of galaxies can achieve a much tighter correlation with 
fc > 1.0. We discuss this point later. 

4. Discussion and Conclusions 

4-1. The Plane is Fundamental 

We have confirmed the conclusion by Koda et al. 
(2000a) that the galactic radius is a fundamental param- 
eter which should be considered to form the best-fitting 
plane in three-parameter space. As one can see from ta- 
ble 1, except for the cases in which both the optically 
measured velocities (V2.2 or {V)) and the isophotal radii 
are used together, the /? values are always around —2.5, 
while the values of a are always around —5.0. Thus, 
approximately, M = — 5.01ogF — 2.51ogi? -|- 7. Since 
M ~ — 2.51ogL, wc have 

L = cV^R. 

Here, c is a constant related to 7, which greatly varies for 
different color bands or different expression of luminosi- 
ties; V^R is expected from the virial theorem. The total 
mass, including the contributions from both the disk and 
the dark matter, should satisfy M{<r) = [V^otl'') '"l/^i so 
that the rotational velocity, Vrot (r) , at radius r is related 
to the total mass, Ai, inside the radius, where, G is the 
gravitational constant. Therefore, our results imply that 
the total luminosity is proportional, on average, to the to- 
tal mass inside the galactic radius. We conclude from an 
observational point of view that the plane represented by 
L ^ V^ R is fundamental for spiral galaxies. 

The mass of spiral galaxies is mostly in the form of 
dark matter, probably in the halo, while the luminosity 
is given mostly by the number of stars in the galaxy disk 
(see figure 5). Therefore, the fundamental plane derived 
by us implies a tighter relation between the bright mass 
in the disk and the dark matter in the halo; at least they 
have similar distribution profiles. This is coincident with 
the conclusion derived from the universal rotation curve 
by Salucci et al. (1993) and Persic et al. (1996) on the 
tight coupling between the dark and the luminous mat- 
ter in spiral galaxies. In general, the luminosity can be 
expressed as L = (TOd/T)A1, where m^ is the fraction of 
disk mass in the total mass, and T is the mass-to-light 
ratio. Our result seems to favor the constant m^ and T 
values for spiral galaxies, which are often assumed in the- 
oretical studies (see Mo et al. 1998; Shen et al. 2001). In 
addition, wc noted that the best-fitting planes (the k val- 
ues) for spiral galaxies in our selected samples vary only 

■^ Courteau (1996, 1997) used the Hubble constant of Hq = 70 
km s~'^Mpc~'^ to calculate the absolute size. This is fine for our 
discussion without further correction. 



slightly in different wavebands, which may result from the 
formation history of the galaxies, and can be used to con- 
strain theories of galaxy formation. 

4-. 2. Tully-Fisher Residuals and Galactic Radius 

Note that the isophotal radii are more effective than 
the disk scale-length in reducing the scatter. As can 
be seen from table 1, the photometry- limited observa- 
tional parameters, such as /J/23.5 from Han (1992) and 
Palunas and Williams (2000), i?H23 from Courteau (1996), 
and Rb25 for the Hubble calibration galaxies, can work 
better. However, when the optically measured veloc- 
ity is used together with the isophotal radii, the best- 
fitting planes are at k > \. Two independent data-sets of 
Palunas and WiUiams (2000) and Courteau (1996,1997) 
show qualitatively the same result. Although we noticed 
the fact that the two observational parameters, log(V2.2) 
and log(i?23.5), are measured at distinctly different radii, 
this may not be the reason for the flat variation of x^ 
when k > 1. 

The fundamental plane involving the linear galactic size 
can reduce the Tully-Fisher residuals by about 50%. The 
unified plane L ex [VRY-^ that Koda et al. (2000a) found 
from the /-band data-set can also work much better than 
the Tully-Fisher relation, though it may not be the funda- 
mental plane. Our results (in our table 1) are consistent 
with case c of theoretical simulations of Shen et al. (2001), 
in which they obtained the smallest scattering (see their 
table 1). 

4.3. Other Issues 

The fundamental plane docs not help to improve the 
distance estimates of galaxies using the luminosity, be- 
cause the galactic size here is a quantity having a dis- 
tance dependence. In fact, during the above searches for 
the plane, the effect as well as uncertainty of the galactic 
distance was diminished for the fundamental plane, com- 
pared with that in the Tully-Fisher relation. 

Based on the Tully-Fisher relation {L ex V^) and 
Freeman's law (R oc i"'^), the fundamental plane de- 
scribed by equation (3) requires a and /3 to satisfy the 
relation a + 2(3 = —10. If a = /3 (see discussion by Zou, 
Han 2000), one may find that a = /3 = —10/3 and ex- 
actly naturally becomes L oc (VR)^^^, identical to that 
obtained by Koda et al. (2000a). Whenever a j^ (3, the 
possible values of a and (3 are in the ranges (0, —10) and 
(0, -5), respectively. The relation of Willick (1999) ob- 
tained from the /?-band, vtf oc L^-^^/g-^^, if converted to 
our manner, is L ~ V'^^'^R^^"^. As can be seen from the 
search results in table 1, these values slightly depend on 
the color band and the sample selection. It is interesting 
to note that these best-fitting a and (3 values are always 
around the central values in their possible ranges. Both of 
the planes derived by Koda et al. (2000a) and by us can 
recover the Tully-Fisher relation and the Freeman's law, 
which are just projections onto 2-D of the fundamental 
plane in 3-D. 

In the above analyses, the bulge and disk parts of the 
spirals were not separately considered. However, for late- 
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Fig. 5. Bulge to the disk luminosity ratio. Most late-type 
spiral galaxies do not have a very bright bulge. (Data from 
Palunas, Williams 2000). 

type spiral galaxies, the bulges often have rather smaller 
luminosities than the disks (figure 5). In our analysis, 
most galaxies are of the late-type, and thus the effect of 
the bulge component on our results is fairly small. The 
disk contribution to the luminosity is dominant. 

4.4- Conclusions 

Our search results reveal that there exists a fundamen- 
tal plane for spiral galaxies, which can be expressed by 
L ~ V^R. The galactic physical size should be involved to 
form the fundamental plane in three-dimensional space 
of logi , logR, and logT^. The size should be prefer- 
ably the linear radius at a given isophotal limit if the 
velocity is taken from the width of H i gas. Otherwise 
the scale length of optical disk can be used as R if the 
rotation velocity is measured from the optical spectrum. 
This fundamental plane can reduce the residual of the 
TuUy-Fisher relation by an amount of about 50%, imply- 
ing that only the other 50% is attributed by the measure- 
ment uncertainties. The fundamental plane exists in all 
optical bands. Such a fundamental plane is probably re- 
lated to the mass and mass distribution of spiral galaxies, 
and should be used to test theoretical work concerning 
galaxy formation. 
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